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Abstract 

We derive asymptotic formulae for central polynomial coefficients, a generalization of binomial coef- 
ficients, using the distribution of the sum of independent uniform random variables and the CLT. 
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For a nonnegative integer ra, Stirling's formula 
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where e is Euler's number, yields an approximation of the central binomial coefficient ( r ^ l / 2 ) (assuming that 
ra is even) using (™) = fc!( ^l fc) , as 
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where we write a m ~ b m as a short-hand for lin^^oo P^ = 1. In the current note, we derive an asymptotic 

formula for the central multinomial triangle or, polynomial, coefficient (cf. [T], [2]), where multinomial 

triangles are a generalization of binomial triangles, where entries in row k are defined as coefficients of the 

{S| polynomial (1 + x + x 2 + . • • + x l ) k for I > 0. Our derivation is not based upon asymptotics of factorials but 

CN| upon the limiting distribution of the sum of discrete uniform random variables {J 

Polynomial coefficients are important, for example, because they denote the number of integer compo- 
sitions of the nonnegative integer n with parts in the set A := {a, a + 1, . . . , 6}, where a, 6, < a < 6, are 
nonnegative integers (cf. [3 ), i.e. the number of possibilities to write n as a sum of integers pi, . . . ,£>&, with 
Pi 6 A, for i = 1, . . . , k. 



2 Multinomial triangles 



In generalization to binomial triangles, (I + l)-nomial triangles, I > 0, are defined in the following way. 
Starting with a 1 in row zero, construct an entry in row k, k > 1, by adding the overlying (/ + 1) entries in 
row (k — 1) (some of these entries are taken as zero if not defined); thereby, row k has (kl + 1) entries. For 
example, the monomial (I = 0), binomial (I = 1), trinomial (I = 2) and quadrinomial triangles (I = 3) start 
as follows, 
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throughout, we assume that m or I are even. If this is not the case, replace respective quantities, e.g. ^r, with their floor, 

ml I 



2 



In the (Z + l)-nomial triangle, entry n, < n < kl, in row /c, which we denote by (V , has the following 
interpretation. It is the coefficient of x n in the expansion of 



(l + x + x 2 + ...+x l ) k = jr( k ) x n . (2.1) 

n =o ^ 1+1 



In recent work, Eger [3] has shown that ( n ) /+1 denotes the number of integer compositions of the nonnegative 
integer n with k parts pi, . . . ,£>&, each from the set {0, 1, . . . , Z}, and allows the following representation, 
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where 7 , is a multinomial coefficient, defined as , , , , , for nonnegative integers ko,...,ki. We 



can easily verify representation (2.2) by noting that for real numbers xq, . . . ,£/, it holds that (cf. [5]) 



(x + xi + ... + x0 fc = Yl 
Then setting a^ = x % for z = 0, . . . , Z, 
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(l + a: + ^ + ...+2: ( ) = > , . U u - fco+ - +( - fci , (2.3) 
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and comparing coefficients of the right-hand sides of (2.1) and (2.3) leads to (2.2). 



3 Generalized Stirling's formula 

Our strategy for deriving approximation formulae for polynomial coefficients is as follows. First, we determine 
the asymptotic distribution of the sum of discrete uniform variables, which we easily find to be a normal 
distribution by the Central Limit Theorem (CLT). Then we determine the exact distribution, which turns 
out to yield the normalized polynomial coefficients (™), . This distribution has recently been found by 
Caiado and Rathie [1] using characteristic functions, residue theory and Cauchy's theorem, but we give an 
elementary two-line proof that almost requires no effort at all. By relating the asymptotic distribution to 
the exact distribution, we obtain a multinomial analogue of Stirling's approximation to the central binomial 
coefficient. 

3.1 Distribution of the sum of discrete uniform variables 

Let m be a positive integer and let I be a nonnegative integer. Let Sj, j = 1, . . . ,m, be identically and 
independently distributed random draws from the discrete uniform distribution on the set {0, . . . , Z}, and let 
g(m) i^ th e i r sum , 
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Obviously, the mean and variance of S^ are given by 

ml (1 A- 1 \ 2 — 1 
Mm>J := E[S< ro >] =mE[S j ] = ™, < ( := Var^™)] = mVarfo-] = ™^^ , 



by standard moments of the uniform distribution. Moreover, by the CLT 
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so that for large m, 5^ is approximately normally distributed with mean ^ and variance m - + 1 ^ — . 

However, we can also exactly determine the probability that S^ 171 ' takes on the integer value n, for 
< n < ml. For this purpose, interpret all Sj, j = 1, . . . ,ra, as independently and identically multinomially 
distributed with probabilities po — • • • — Pi — jh[ of 'events' to Z occurring, i.e. each 5j is vector- valued, 
Sj = (Aq, . . . , A/), with P[Sj = (ao, . • . , a/)] = ^rj for ao + . . . + a\ = 1, a& > 0, where Ak denotes the 
number of times event k occurs, for k = 0, . . . , I. Hence, the sum S^ is multinomially distributed with 
probability mass function 

P[S^ = (a ,...,a l )}=P[A =a ,...,A l =a l }=( ™ )(^-f, 
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where ao -\- . . . -\- ai = m. Then, if S^ = (ao, . . . , aj) for nonnegative integers ao, . . . , aj with ao + . . . + a^ = m, 
g(™) 'represents' the integer • ao + ... + /• a/. Thus, for n such that < n < raZ, we obviously have, 
reinterpreting S^ as scalar-valued, 

P[s<m '^ ^S_ P[S, " >M "° °' H ^ r Jx- L".J-(<ti)"(-)« 

0-ao + ...+^-az=n 0-ao + ...+£-az=n 

3.2 Approximation Formulae 

From our above derivations, we know that the normalized random variable S^ converges in distribution 
to a normally distributed random variable, which precisely means that the distribution function of S^ 
converges pointwise to the distribution function of the corresponding normally distributed random variable. 
Let < n < ml and consider 

P[S (m) = n] = P[S (m) < n] - P[S (m) < n - 1]. 

By convergence of distribution functions, we hence know that, for large m, 

P[S (m) = n] w P[n - 1 < X < n] « P[n - 1/2 < X < n + 1/2], (3.1) 

where X has a normal distribution with the same mean and variance as S^ m \ and applying a continuity 



correction (cf. [4]). Rewriting the right hand side of (3.1) in terms of the standard normal variable Z 
~ MrrM , we thus obtain for n = a m i , 

P[S^ = ^] « P[--J- < Z < +-J-] = $ (^-) - $ (--1 

where $ is the standard normal cumulative distribution function. Geometrically, <l> f ^r — ) — ^ ( — 2^ — ) 



denotes the area under the normal curve <p(z) = l/>/27rexp(— z 2 /2) from — -^ — to +2^ — • ^ or ^ ar S e m > 
this area can be approximated by the area of the rectangle with width -^— and height 0(0), that is, 

P[S(m) = —] « X 
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Using the expression for P[S^ m ' = n] from above, we hence have 
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For 1 = 1, Pascal's case, I = 2, I = 3, and Z = 4, we therefore have the approximations 
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Figure 1: Distributions P[S^ m ^ = n) for m = 5, 10, 20 for I = 4 fixed; and normal approximations. 



In Figure 3.2, we show for I = 4 the distributions P[S^ = n] for m = 5,10,20 and their respective 
normal approximations. 

These plots suggest to more generally consider a pointwise approximation of P[5^ m ) = n], for < n < ml, 
by the value f(n) of the density of the corresponding normally distributed random variable X, leading to 
the following approximation formulae for polynomial coefficients (^), , 
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How well are these 'normal approximations' a m?rij / of a m ^ n j := (^) /+1 for various values of n for finite ml In 
Figure 3.2, we depict the relative error I am > n > l ~ a " m > n > 1 I for all values of n between and ml for m = 100, 200 

I Q"m,n,l i 

and / = 3,4. We see that the relative error becomes very low for n 'around' ^ and increases symmetrically 
as n approaches the boundary values and ml. In the same figure, we show the decrease of the relative error 
for the central polynomial coefficient as a function of m. 
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Figure 2: Left: Relative error of normal approximation of (™) for m = 100, 200 and I = 3, 4 as a function 
of n. Right: Relative error for central polynomial coefficient as a function of m for I = 4. 



